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Abstract. We introduce a class of Markov stochastic processes called m- 
polynomial, for which the calculation of (mixed) moments up to order m only 
requires the computation of matrix exponentials. This class contains affinc 
processes, Feller processes with quadratic squared diffusion coefficient, as well 
as Levy-driven SDEs with affine vector fields. Thus, many popular models such 
as the classical Black-Scholes, exponential Levy or affine models are covered 
by this setting. The applications range from statistical GMM estimation to 
option pricing. For instance, the efficient and easy computation of moments 
can successfully be used for variance reduction techniques in Monte Carlo 
simulations. 



Pricing and hedging of contingent claims is the crucial computation done within 
every model in mathematical Finance. For European type claims this amounts to 
the computation of the expected value of a functional of the (discounted) price 
process under some martingale measure. Hedging portfolios are then constructed 
via appropriate derivatives of those expected values with respect to model param- 
eters or to the prices, so called "Greeks". Let us denote the (discounted) price 
process at time T, a vector in W l , by Xp. We can roughly distinguish three cases 
of complexity for the mentioned computations: 

(1) The probability distribution of Xp is known analytically. 

(2) The characteristic function of X? is known analytically. 

(3) The local characteristics of Xt are known analytically. 

In the first numerical quadrature algorithm is sufficient for the efficient 

computation of the contingent claim's price E[<P{Xt)\, where 4> denotes some payoff 
function. In the second case, variants of Plancherel's theorem are applied in order 
to evaluate the price E[<P(Xt)], for instance, 



where <j> denotes the Fourier transform of the function <f>. Remark that often mod- 
ifications of the original payoff function are used to make the Fourier methodology 
applicable. This is numerically efficient, even though its implementation (in par- 
ticular the complex integration) can take some time (see, e.g., [3]). Also there are 
different levels of what it means to "know analytically" the characteristic function 
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of Xt- If one, e.g. in affine models, has to solve a high-dimensional Riccati equation 
for each u £ M. n to calculate the characteristic function u i— > E[exp(i(u, Xt))}, the 
"analytic knowledge" means at least some pre-calculations which have to be per- 
formed efficiently, too. In other words, methods applying the characteristic function 
are efficient if the amount of pre-calculations is limited. 

The third case is characterized by the use of Monte-Carlo simulation methods: 
one samples from the (unknown) distribution of Xt by generating, for instance 
through Eulcr schemes, approximate distributions for Xt- From those approximate 
distributions, which should be easy to simulate, one draws a sufficient amount of 
samples. This two-step procedure is very robust, but takes a considerable amount 
of time. 

In this article we would like to add a fourth case, which in the previous order 
would correspond to case 2~. We can describe a class of processes, called "poly- 
nomial processes" , where it is easy and efficient to compute moments of all orders 
of the random variable Xt, even though neither the probability distribution nor 
its characteristic function need to be known. We shall analyze this class and show 
that exponential Levy processes, affine processes or Jacobi-type processes belong 
to it. The method is best explained by an example: consider a stochastic volatility 
model of SV J J- type [5], i.e. both the logarithmic (discounted) price processes and 
the stochastic volatility can jump. Such models can be described by stochastic 
differential equations of the type 

dX t = (aX t + bV t + c)dt + \fV t dBt,i + dZ tA , 
dV t = (aV t + f3)dt + Vv t dB t:2 + dZ t>2 , 

where (-61,-62) are possibly correlated Brownian motions, and (Z l7 Z 2 ) are Levy 
processes, independent of (-61,-62), where the second component Zi is a subordi- 
nated For such models there is no easy-to-implement (explicit) formula for the 
characteristic function, even though they arc affine models. Assuming now appro- 
priate moment conditions for the jump measures, the Markov process {X, V) has 
the remarkable property that the expected value of any polynomial of the process 
(X, V) is again a polynomial in X and Vq. The coefficients of this polynomial can 
be calculated efficiently by exponentiating a matrix, which can be easily deduced 
from the generator. In other words, there is a dense subset of claims in the set of 
"all" claims, where the prices and hedging ratios are known explicitly (up to matrix 
exponentials). This explicit knowledge allows to compute prices of general claims 
by variance reduction techniques, which is more efficient (and easier to implement) 
than Fourier methods with "unknown" characteristic function. 

The article is devoted to a study of polynomial processes in some depths together 
with their most obvious applications to mathematical Finance. As the most striking 
application we see the possibility to reduce the variance of Monte Carlo evaluations 
within a polynomial model by approximating the claim through elements of the 
dense subset of already priced claims. The remainder of our article is organized as 
follows: in Section [2] we introduce formally the class of m-polynomial processes and 
draw several basic conclusions. In|3]we analyze the conditions on the characteristics 
of a Feller semimartingale to be m-polynomial. Section [4] deals with examples from 
the class of m-polynomial processes and Section [5] with applications to pricing and 
hedging in mathematical Finance. 
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2. Polynomial Processes 

We consider a time-homogeneous Markov process X := (Xf) t>0 xeS with state 
space S C E™, a closed subset of R™, and semigroup (Pt)t>o defined by 

PJ(x) := E[/(Xf )] - / /(OPt(z, de), x e S, 
Js 

and acting on all functions / : S — > E which are integrable with respect to the 
family of Markov kernels p t (x,-). 

Notation 2.1. Le£ Pol< m (S) denote the finite dimensional vector space of poly- 
nomials up to degree m > on S 1 , z.e £/ie restriction of polynomials on E™ to S, 
defined by 

»ki i e S, Ok e K > , (2.1) 

[|k|=0 J 

where we use the multi-index notation k = (fci, . . . , G Nq, |k| = fci + • • • + k n 
and Pol< m (S') is endowed with some norm \\ ■ ||p i <m a,nd its 

dimension is denoted by N < oo. Moreover, Pol m (S') corresponds to the vector 
space of polynomials which are precisely of degree m. 

Definition 2.2. We call an S-valued time-homogeneous Markov process 
m-polynomial if 

PJ(x) e Pol< m (5) 

holds true for all f € Pol< m (5') and t > 0. If X is m-polynomial for all m > 0, 
then it is called polynomial. 

Theorem 2.3. Let X be a time-homogeneous Markov process with state space S 
and semigroup (Pt), pointwise continuous at t = 0, i.e. (h Ptf(x) is continuous 
at t = for all x e S and f : S — > E where P t f exists as finitely valued function 
on S. Then the following assertions are equivalent: 

(i) X is m-polynomial for some m > 0. 

(ii) There exists a linear map A on Pol< m (S l ), such that (Pt) restricted to 
Pol< m (S') can be written as 

^t|Pol< m (S) = e 

for all t > 0. 

(iii) The infinitesimal generator A is well defined on Pol< m (S') and maps 
Pol< m (S) to itself. 

(iv) The Kolmogorov backward equation for an initial value /(•, 0) G Pol< m (5) 

d -^ = Af( X ,t) 

has a real analytic solution for all times t € E. In particular, f(-,t) € 
Pol< m (S). 

Proof. We start by proving (i) ^> (ii). Let £(Pol< m (S)) denote the space of all 
linear maps on Pol< m (5). By the semigroup property, 

P ( .)|poi< m (S):R+-APol< m (S)) (2.2) 
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satisfies the Cauchy functional equation 

frvy all -f o ~> n 

(2.3) 



P t+S = p t p s for all t, s > 



Po = I 



Since Pol< m (5) is finite dimensional, continuity of (2.2 1 at t — already implies 
that there exists an A S £(Pol< m (5)) such that -Pt|p i <m (s) = e*" 4 - 
Next, we show (ii) => (iii). By the definition of the infinitesimal generator, we have 
for all / G Pol< m (5) 

Ptf-f e tA f - f 

Af = lim ^ = lim J ^- = Af, 

t^o t t^o t 

which is obviously well defined with respect to || • ||p i <m and in Pol< m (5). 
Assertion (iv) follows from (iii) since A maps Pol< m (S f ) to itself. The Kolmogorov 
backward equation can therefore be understood as a linear ODE in the classical 
sense, whose solution is thus given by e tA f(-, 0). 

At last, we prove that (iv) implies (i). For any initial value / in an appropriate 
Banach space, P t f is the unique solution of the Kolmogorov equation. On Pol< m (<S') 
it must therefore be equal to e tA f and is hence a polynomial of degree smaller than 
or equal to m thereon. Thus, X is ro-polynomial. □ 



Remarks 2.4. (i) There is no need in Theorem \2. 3jf ii) to restrict the time param- 
eter t to R + . Actually, for t G R, (e tA ) extends to a group generated by A. 

(ii) The assumption of pointwise continuity of Ptf at t = can be replaced by differ- 
ent conditions. For example, stochastic continuity of X together with the existence 
of moments up to order m + e, for some e > 0, is another sufficient hypothesis for 
Theorem \2.3l 

(iii) It is certainly possible to consider other finite dimensional functional spaces 
instead of Pol< m (5'). For example, trigonometric polynomials are left invariant 
by Markov processes of the form Xf — x + L t , where L t is a Levy process. We 
have been choosing polynomials since we wanted to enlarge the class of affine pro- 
cesses. This statement is not true literally, since affine processes which do not admit 
first moments, can never be polynomial, however, if a regular affine process admits 
moments of order m it is m-polynomial. Therefore properties of m-polynomial pro- 
cesses apply to affine processes which admit moments up to order m (see Example 



4.1) 



Corollary 2.5. Let X be an m-polynomial process with semigroup (Pt), continuous 
at t = and let f G Pol< m (5') be fixed. Then there exists a unique function 
Q : K x S — > M., being real analytic in time and Q(t, •) € Pol< m (S') for all t 6 R, 
such that 



Q(0,x) = f(x) (2.4) 



and 



Q(t — s, X s ) is a martingale for s > 0. (2-5) 
Proof. By the martingale property of Q(t — s, X s ), we have for all s > 

Q(t,x)=E[Q(t-s,X;)]. 
Thus, by setting s — t we obtain 

Q(t,x) = E[Q(0,Xf)] = E[/(Xf)] = Ptf(x), 
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since Q(0, Xf) — f{Xf). As X is m-polynomial, Q(t,-) e Pol< m (S) and by the 



previous Theorem it is a real analytic function. Clearly, P t f(x) satisfies (2.4) 



and p.5\ . □ 

We recall the following definition (see e.g. [IT], [H]). 

Definition 2.6. A real valued function R : [0, oo) x S — > M is said to be time- 
space harmonic for X if (R(s, X S )) S > Q is a martingale. If moreover, R(s,x) is a 
polynomial in s and x that is time-space harmonic for X, it is called a time-space 
harmonic polynomial. 



For t = 0, Q{— s,X s ) as defined in Corollary 2.5 can be considered as a time- 
space harmonic function for the m-polynomial process X. If A(Pol< m (S)) C 
Pol<( m _ 1 )(5'), then it is a time-space harmonic polynomial, which is proved in 
the Corollary below (see also [T] for S = K) . 

Corollary 2.7. If A(Pol< m (S)) C Pol<( m _i)(5), then Q(—s,x) as defined in 
Corollary \2.5\ is a time space harmonic polynomial. 

Proof. For some fixed / G Pol< m (S) we have 

Q(-s,x) = P_J(x) = e~ sA f(x) = V ^-A n f{x) = V { -=^-A n f{x), 

TV. TV. 

n=0 n=0 

which is obviously a polynomial in s and x and a martingale due to Corollary |2.5| 

□ 

3. Polynomial Feller semimartingales 

For the class of Feller processes, we want to find sufficient conditions for m- 
polynomial processes in terms of the infinitesimal generator. We consider therefore a 
Feller semigroup (P t ) on S. If C^°(5) is contained in the domain of the infinitesimal 
generator A, then it is well known that there exist real- valued functions aki,bk,c 
and a kernel K(x, d£) on S x 6(M"\{0}) such that for u G C^(S), A is given by 

Au(x) = I V dki (x) f U ^ + J^bk (x) - c(x)u(x) 



k,l — l k—1 



J 

J'R 



^ / du(x) 



(3.1) 



"\{o> V fc=l 



where a(x) = (dfei(aj))/ !) i=i,... 1 n is a symmetric positive semidefinite matrix, 6(a;) G 
K™, c is non-negative, if (x, •) is a Radon measure onSx S(R™\{0}) and x : - * 
K n some bounded continuous truncation function with x(0 = £ in a neighborhood 
of (see for example [TT] or [TB]). Clearly, the above parameters have to satisfy 
certain admissibility conditions guaranteeing the existence of the process in S. 
As the definition of m-polynomial processes requires the existence of moments up 
to order m and thus the finiteness of X t a.s., we henceforth assume X to be con- 
servative, i.e. c = 0. 

If X is additionally a semimartingale, which is automatically the case if S — W l 
(see Theorem 7.16 in [4]), then its characteristics (B,C, v) associated with the 
truncation function x(£) are given by 

B t = ( b(X s )ds, C t = [ a(X s )ds, v{dt,d£) = K(X t ,d£)dt. 
Jo Jo 
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We shall refer to (b,a,K) as differential characteristics of X (see |14j). 

In order to specify the form of a, b and K such that A generates an m-polynomial 

process, we start with the following Lemma. 

Lemma 3.1. Let Xf be a conservative Feller process on S. If for every f G 
Pol< m (S) , there exists a polynomial g G Pol< m (S f ) and some to > such that the 
process 

C{ = f{X't) - f(x) - [ g(Xf)ds (3.2) 
Jo 

is a well defined martingale for t G [0,to]j then X is m-polynomial. 
Proof. Due to the martingale property, we obtain for t < t 



P t f(x) - f{x) 



P s g{x)ds = E 



f(Xf)-f(x) 



g{X:)ds 



0. 



Let us denote by Sg = SU{6}, (5 ^ S) the Alexandrov one-point compactification of 
S. By the Feller property of (P t ), PJ G C(S S ) for aU t > Q and || • = sup xeSs \ ■ \ 
is well defined. Furthermore, for t < to 



(Ptf -f)~9 



1 f* 

7 / (P s g-g)ds 
1 Jo 



< 



\Ps9 



ids, 



which goes to as t goes to since by the Feller property lim s ^ ||-Ps<7 — <?||oo = 0. 
Hence, with respect to || ■ ||oo, Af is well defined. Thus, it is also well defined with 
respect to |j • ||p i< m (s) an d Af — g G Pol< m (S'). As P t is continuous at t = (by 
the Feller property), we can apply Theorem |2.3| By assertion (iii) it then follows 
that X is m-polynomial. □ 



Remark 3.2. The operator Q satisfying g = Qf in equation (3.2) for f G D(G) is 
usually called the extended infinitesimal generator. 

In the next Lemma we establish a maximal inequality for semimartingales which 
lie - as we shall show in Thcorcm |3.6| below - in the class of m-polynomial processes. 
A similar statement and proof for the case of Levy driven SDEs can be found in [T2"] . 
see also [TS]. Before stating the Lemma, let us define the following two conditions 
on the kernel K(x,d£): 



Condition A. The kernel K(x,d£) is of the form 

K{X t ,dt):=tioo(d§ + J2 x s,i»io(d£)+ J2 WjWj^). 
ie/ (i,j)e,J 

where all foj are Levy measures on R n with 



€11 >i 



(3.3) 



-} 



The index sets I and J are defined by 

I = {1 <i < n\S t C 
and 

J = i < j\Si x Sj C R 2 + or S t X Sj C 

where Si stands for the projection on the i th component. 



-}, 
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Condition B. The kernel K(x,d(;) satisfies 

K(X t ,d£) :=g?*ti(dt), 

where for each x € S, g x fi denotes the pushforward of the measure fi under the map 
g x . Moreover, g x {y) — g{x,y) is affne in x, i.e. 

5:5x1^1", (x,y)»H(y)x + h(y), 

with H : M. d — > M. nxn and h : M. d — > M. n some measurable functions. Finally, fi is a 
Levy measure on M. d integrating 

{\\H(y)\\ k + \\h(y)\\ k ) nidy) for all 1 < k < m. (3.4) 



Remark 3.3. The integrability condition (3.4 1 (in particular for k = 1 ) can be 
weakened for specific choices of H and h (see Example 4-3). 

Lemma 3.4. Let m > 2 and Xf be a semimartingale with the following differential 
characteristics (b,a,K) associated with the "truncation function" x(0 = Q 

n 

b t = b + J2 x t,ifc, b, I3i e R", 



a t = a + ^X tii a i0 + .\,.,.Y».,n,,. a, a„ S R nxn , 

i—l i<j 

and the kernel K satisfying either Condition A or Condition B. Then, there exists 
a constant C m such that for every t £ [0, 1] , 

ft 



E 


su P rar 


<c m 






_s<t 







jE[\\X„\r]daY (3.5) 



Proof. Due to the choice of the "truncation function" x(0 = the canonical 
representation of X is given by 

t(n x (d8,dt)-v(da,d£)) 



Xf = x + B t + X c t 
= x + B t +X. 



•\{o} 



t t ' 



where fi denotes the random measure associated with the jumps of X, X c the 
continuous local martingale part and X d the purely discontinuous martingale part. 
Using the inequality 

\\x + y\r<2 m - 1 (\\x\r+\\y\n (3.6) 

for m > 1, we have 



su P |ix s a; ir<4 

s<t 



m— 1 



sup \\B„ 



-sup|ix s c ir+su P ||x s c 

s<t s<t 



Thus, in order to prove (3.51, it suffices to find estimates for B, X c and X . For 
the sake of notational simplicity we only consider the one dimensional case in the 
following computations. (We also write | • | instead of || • ||.) Moreover, the constant 



^The integrability assumptions on the Levy measures constituting the compensator u allow 
the choice of this "truncation function" . Here, /||£||>i IICIlMij {dQ < 00 m case of Condition A and 
JK d \{o}(ll^(^)H IIM^IDM 6 '!/) < oo in case of Condition B are sufficient. 
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C m in the inequalities below may vary from line to line and t is always assumed to 
be in [0, 1]. We start with B: 



sup|B s | — sup 

s<t s<t 



{b + XuPi) du 



< sup ( [ \b + X u (ii\du 
s<t \Jo 



<n \b+ x u p x \dv)j =(j^\b+x u 



Pi\ t{u<t}{u)du 



< / \b + X u p l \ m {l {u < t} {u)) m du<C n 



+ \X u \ m \p 1 \ m )du, 



where the penultimate inequality is a consequence of Jensen's inequality and the 
last one follows from (3.6). Fubini's theorem finally leads to 



E 



sup \B S 

s<t 



<C m I (l+E[\X s \ m ])ds. 
Jo 



In the case of X c , an application of the Burkholder-Davis-Gundy inequality yields 



E 



sup\X c s 

s<t 



< C m E 



[X\X% 



By our assumption, we have [X c ,X c ] t = J Q (a + X s a w +I s 2 a n ) ds > 0. Hence, 
we obtain 



(a + X s aiQ + Xg an) ds < I J \a + X s aio + X 2 S an | ds 

< [ \a + X s a w + X^an\^ (l {s < t} (s))^ds 
Jo 

<C m [ (|a|f +\X.\ 3 t\ai Q \¥ + \X.\ m \a 1 i\ 3 *)ds, 



where the two last expressions follow again from (3.6 1 and Jensen's inequality since 
we supposed m > 2. Thus, we have 



E 



sup \X C B 

s<t 



<C m [ (1 + E[|X S |?] +E[\X s \ m ])ds. 
Jo 



We finally turn to the purely discontinuous martingale part X d . A further appli- 
cation of the Burkholder-Davis-Gundy inequality yields 



E 



sup \X° S 

s<t 



<C m E 



[X d ,X d ] ] 



where [X d ,X d ] t = £ < s < 4 \ AX ^\ 2 ( see for example Q3], 1.4.51). Since B and X c 
are both continuous, we therefore have 



Y t := [X d ,X\ = \ AX s\ 2 = f I 

Jo Jm. 



0<s<t 



'0 Jm.\{0} 



\tf» x (ds,dO. 
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Following the approach of [12], let us define T n = inf{< | Y t > n}. Then, 

4„= E (n_+AY.)*-(y.js 

s<*AT„ 
ptAT n p 

= / ((r._ + K| 2 )*-(Y._)*) 

^0 JR\{0} 

which is due to the fact that Y is purely discontinuous, non-decreasing and AY S 
\AX S \ 2 . Furthermore, since v(ds,d£) is the predictable compensator of fi x 









f(AT„ /• 


E 


Y 2 

tAT n 


= E 


Jo Jb 



((y._ + iei a )*-(y._)*)K^.dC) 

The inequalities (see [12] ) 

iy + z )P-yP< 2 P- 1 (yP- 1 z + Z P ), 
for a;, y, z > 0, e > and p > 1, then yield 



(3.7) 

(3.8) 
(3.9) 



(n_ + - < [ ey?_ \e + —\e + ier 

V £ 2 1 



as well as 



((y s _ + iei 2 )^ - ^ < 2^ ( ey.i lei 2 + ^Tiei 2 i^i T + ien*.l 



1 



for X s > 0. In the first case, inequality ( 3 . 9 1 is applied to (Y s _)~2~ 1 , i.e. x = 1 
and in the second case to (Y S _) ! %~~ 1 X S . If the kernel satisfies Condition A, (3.7) 
can therefore be dominated by 

1 











E 


Y 2 

tAT n 


< 2 i s i - 1 E 


TV ( 

JO JR\{0} V 



tAT„ 



JR\{0} 
tAT„ 



JR\{0} 



^n! lei 2 + ^T3iiei 2 l^l f + KH*.I ) MdOds 

S 2 1 I 



£ 2 1 ' 



< 2^- L £ 



K\{0} 



ie| 2 (Moo(^)+Mio(dO + Mn(^)) E 



bt A K 



tAT„ 



+ 2^(tAr„)/ (^rl£| 2 + 

JRUO} V e 2 



+ 2^ 
The inccjuality 



R\{0} 

tAT„ 

(i 

*AT„ 
JR\{0} 



ier moo(^) 



^=rl£| 2 E [|X s |f ] + \t\ m E[\X s \]\ n w (d$ds 



E 



tAT„ 



r 2 ds 



< E 



n 2 A Y 4 



tAT„ 
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follows by the fact that Y s _ is non-decreasing and Y s < n for s <T n . Observe that 
due to the moment assumptions on the Levy measures, the integrals with respect 
to /Ltoo, Mio and /Un are well defined. Choosing for example 



e = 



R\{0} 



|C| 2 (Moo(^)+Mio(^) + Mn(dO) 



leads to 



2 













y ? 
tAT„ 


< E 




- l E 
2 





tAT„ 



(1 + E[|X S |]+E[|X S | 2 ] +E[|X S |^] + E[|J*q m ])ds. 
Since Y t ^ T ^ < Y t * Tn i , we can apply the monotone convergence theorem to obtain 
< C m I (l + E[\X a \]+E[\X s \ 2 ] + E[LY 5 |t] + E[\X s \ m ])ds. 









E 




<c m f 






Jo 



In the case of Condition B, (3.7) becomes 









/>t/\T n r 


E 


y *2 
tAT„ 


= E 


JO JiS 



((Y s _ + \H(z)X s + h(z)\ 2 )^ - (Y s _) t ) M (dz)ds 



Again, by (3.8 1 and ( |3.9| , the integrand can be dominated by 

(F s _ + \H(z)X s + h(z)\ 2 )^ - (Y s _) ? < 



2^ 



(y s ? 1 \H{z)X s + h(z)\ 2 + \H{z)X s + h(z)\ m ^) < 



2^(Y s f x ||H(«)|| J |X.| J + y a ! >(z)| 2 +-|i7(z)X s + ^)| m ) < 
2^ (erZdli/WII 2 + |M^)| 2 ) + -4^(\X s r\\H(z)\\ 2 + \h(z)f 

\ £ 2 1 

2 m - 2 |x ! nii?(z)ii m + 2 m - 2 |/i(z)r 

Note, that the integrability assumptions on /i guarantee 

(||ff(z)|| 2 + |^)| 2 ) n[dz) < oo, 

diff(«)ir+i/i(«)r)/i(<te) <oo. 

'\{Q} 

By applying the same arguments as before and by choosing 



(3.10) 



£ = 2-(? +1 > 



(iiffwr+iMz)! 2 )^) 



we obtain 









E 




5: Cm / 






JO 



K d \{0} 



<C m / (l+E[|X s | m ])ds<oo. 
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Combining all these estimates, we finally have for t £ [0, 1] 



E 



sup \X 3 S 

s<t 



< 



C m [\x\ m + / (1+E[\X S \] + E[\X S \ 2 ] + E[\X S \™] +E[\X a \ m ])ds < 



Cm [ N m + J (l+E{\X s \ m })d S 



□ 



Remark 3.5. The proof of the above Lemma in the multidimensional case is com- 
pletely identical. One only has to make use of the equivalence of the p-norms and 
the fact that | £\ d XiXj | f < C||a:|| m . 

Theorem 3.6. Let m > 2 and Xf be a conservative Feller semimartingale on 



S whose infinitesimal generator on C^°(S) is of form (3.1l with c — 0. Assume 
furthermore thatE )|| m ] < oo for allt £ [0, 1] (of course [0, 1] could be replaced 
by [0,e] for any e > 0). Then, X is m-polynomial if its differential characteristics 
(6, a. K) associated with the "truncation function" x{0 = 6 are as i n Lemma 3.4 
with K satisfying either Condition A or B. 



Proof. Let u s the denote the right side of (3.1 1 literally by In view of 



Lemma 3.1 we show that for every / £ Pol< m (5), f £ Pol< m (S) and that 



M( := f(xn - f{x) - / A*f(X s )ds 
Jo 



(3.11) 



is a well defined martingale for t < 1. Indeed, our assumptions on a, b and K, 
imply A# (Pol< m (S)) Q Pol< m (5). This is obvious for the differential operator 
part of A& . For the integral part, which we denote by Af , we have in the case of 
Condition A 



Aff(x) 



>\{0} 



fe=i 



dx k 



E / /(*+o-/(*)-£e 



fc — — I Xifj,i {d€) 



fc=i 



(4j)eJ ^R"\{0} V fc=i 



g/(g) 

dx k . 



XiXj fJ,ij {d^) , 



which lies in Pol< m (5), since /(# + f) - /(a;) - J2k=i £k~B^r is in Pol< m - 2 (5')- 
By the integrability conditions on the Levy measures /Ay the integral is also well 
defined and the truncation function \ can be chosen to be the identity. In the case 
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of Condition B, Af is given by 



Affix) = / (f( x + ^-f( x )-J2i 
[ (f(x + g(x,y))-f(x)-Y,(g(x,y)) 
[ (f{x + H(y)x + h(y)) - f(x) - f^{H (y)x + h{y)) 



df(x) 
dx k 

df{x) 
dx k 



K(x, d£) = 

fi(dy) = 

df(x) 



dx k 



H(dy). 



IR«\{0} y fe=1 

This is clearly in Pol< m (5) and by the assumptions on the Levy measure /i, namely 
\H{y)\\ k + \\h{y)\\ k ) fi(dy) < oo, for aU 1 < k < m, 

R d \{0} 



the integral is also well defined. Hence, only the martingale property of (3.11) 
remains to be shown. By applying Ito's formula it follows that M{ is a local 
martingale. However, since we have A#f £ Pol< m (S') for / 6 Pol< m (S f ), we can 
dominate sup s<t A// for t < 1 with the following expression: 



sup ||m/|| = 



sup 

s<t 



f(X?) ~ f{x) - A*f{XZ)du 



< sup || f(Xf) || + || f{x) ||+ sup 

s<t s<t 



A*f{XZ)du 



< d ( i + sup \\xf\r ) + c 2 (i + \\x\n + 1 c 3 (i + \\xf\n d S 

s<t 



< Ci 1 + sup \\X*\\ m + C 2 (1 + \\x\\ m ) + C 3 1 + sup \\X?\\ 



< c + \\x\r+ S up \\x*\r^j 



with C — C\ + C 2 + C3 for some constants Ci, C2, C3. By Lemma 3.4 and our 
moment assumption on X, we know that E [sup s<t \\X^ || m ] is finite for every t £ 
[0, 1]. This implies that 



E 



sup 

s<t 



\m!\ 



< 00 



for every t < 1. Thus, M{ is a martingale for t £ [0, 1] and Lemma |3 . 1 1 yields the 
assertion. □ 

Remark 3.7. It is important to note that the differential characteristics of X in 
Theorem 3.6 are specified with respect to the "truncation function" = £. While 



a and K do not depend on the choice of x> the characteristic b — b(x) does. So, 
if one chooses another truncation function x instead of then b{x) transforms as 
follows 



Mx(0) = MO 



'\{0> 



m)-OK{x t ,do. 
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Thus, the requirement that a and K are as in Theorem \3.6\ ano 
is an equivalent condition guaranteeing that X is m-polynomial. 



[ (Z-m)K(X t ,do) ePol<!(5) 



Remark 3.8. Note that an operator of form (3.1) with x(0 = £ can satisfy 

A(Pol< m (S)) C Pol< m (5), (3.12) 
even though the conditions of Theorem\3.6\ are not fulfilled. In the case of m = 2. 



in order to satisfy (3.121, K(-,dt;) can be any function as long as 



consider for example the one dimensional process with b as in Theorem 3.6 Then. 

(3.13) 



a(x) 



f ?K(x,dS) ) ePol< 2 (S). 
Jr\{o} J 



If m > 3 and (3.131 is satisfied, then (3.121 implies if(-,<i£) G Pol<3(S f ). 



Remark 3.9. As an assumption in Theorem \3.6\ we required M[\\Xf\\ m ] < oo for 
t G [0, 1]. In the case of Levy SDEs of the form 



X, = x- 




x s ,i a io + x s ^x 

i—1 i<j 



dB, 



(3.14) 



(H(y)X s + h{y)) (u x (ds,dy) - »(dy)ds), 

lO JR d \{Q} 

where B t is standard Brownian motion in R" and fj, (ds,dy) a Poisson random 
measure with compensator ji{dy)ds satisfying (3.4 1, this condition is automatically 
fulfilled if 



a + X Sy iaio + X S iX S j 

i—l i<j 

is Lipschitz continuous on S . This can be shown by proving that the map h x (Z) t 



A x (Z) t = x + 



I (b+Y^Z s ^}jds 




i=l i<j 



dB„ 



{H{y)Z s + h(y)) (fi x {ds, dy) - u(dy)ds), 



is a contraction on the Banach space C pr ([0,t]; L m (f2,]R")) (see [7] for details). The 
proof can be based on the same arguments as those used in Lemma ( |3.4[ ) . A similar 
result is also stated in 1121. 
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In the case of affine characteristics, the existence of the moments of the Levy mea- 
sure in the sense of (3.3) might also be sufficient for K[\\X^\\ m ] to be finite (see 
Lemma 5.3 and 6.5 in 

4. Examples 

In order to apply Theorem |3.6| to the following examples, we assume throughout 
this section m > 2, i.e. the process X will admit moments up to order at least 2. 

Example 4.1 (Affine processes). Every conservative, regular affine process X on 
S = R p + x R n ~P is m-polynomial, j/E[||(Xf)|| m ] < oo for t G [0, 1]. For details on 
affine processes see [5] . 

Proof. On Cg, the generator of an affine process is given by 

1 " ( V \ d 2 u(x) I \ 

Au(x) = - ^ I a k i + ^ Xia i0> kl J q x q- + ( b + ^ x i&' Vtt(i) ) 
k,i=i \ i=i J x k %i \ i=1 J 

s\m V l= i / 

where the parameters have to satisfy certain admissibility conditions (see [5]). By 
Theorem |3.6| and Remark |3.7| it then follows directly that X is m-polynomial. □ 



If a conservative affine process has constant drift and exhibits only state inde- 
pendent jumps, ^4(Pol< m (5)) C Pol<( m _ 1 )(S l ) holds true. In this case Q{~s,x) = 
P- S f{x) is a time space harmonic polynomial by Corollary 2.7 

Example 4.2 (Levy processes). Let L be a Levy process on M." with triplet (b, a, fi) 
satisfying Jjijii^ IICI| m M(^0 < 00 • Then, the Markov process Xf — x + L t is m- 
polynomial. 



Proof. This can immediately be seen from Theorem |3.6| since b, a and fi are simply 
constant functions. □ 

Example 4.3 (Exponential Levy models). Exponential Levy models are of the form 

Xf = xe L \ (4.1) 

where L t is a Levy process on K with triplet (b,a,fj,). X t usually corresponds to 
the price process under a martingale measure. The absence of arbitrage then im- 
poses conditions on fi and b since the discounted price process e~ rt X t must be a 
martingale: 

e y n(dy) < oo and b = r- -a— [ (e y - 1 - yl\ y \<i) fi{dy). 

\y\>l 1 JR 

The infinitesimal generator of ( |4.1| is given by 

. , . ax 2 d 2 u(x) du(x) f ( , ... , . ... „. du(x)\ , , . 

Au(x) = — r^ 1 + rx— — + [u (xe y ) - u(x) - x (e v - 1 — M- (i(dy). 

2 ax z dx . a \ ax J 



In terms of Theorem 3.6 we are in the situation of Condition B with g(x,y) 
H(y)x = (e y — l)x. Note that the integral part of A is well defined for u € Pol< m (I 
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if 



Indeed, by Taylor's theorem 



ay Kdy) < 



(4.2) 



u (xe y ) - u(x) -x(e y - 1) I fi(dy) 



\{e v - \)x\* sup 

\ V \<1 «<8U P | B |<! |(e«-l)x| 



dx 
d 2 u(z) 



< 



dz 2 



Hidy)- 



\C(x)\(l+\( e y -1)D n(dy), 



ll/l>l 



where C(x) g Pol< m (IR). Since \(e v — 1)| 2 < C\y\ 2 for \y\ < 1 and some con- 



stant C and since \(e y - l)| m < 2 m " 1 (e my + 1), w e obtain that (|42j is sufficient 
for the existence of the integral. Condition (3.101 in the proof of Lemma 3.4 is 



then automatically fulfilled. Furthermore, (4.2 1 guarantees also the existence of 



E [| (Xf )\ m ] < oo, whence by Theorem 3.6 exponential models are m-polynomial 
if (4.2) holds true. Note that the Black-Scholes model falls into the realm of this 



example. 

Example 4.4 (Levy driven SDEs). Let L t denote a Levy process on R d with gen- 
erating triplet (b,a,n). Suppose furthermore that Vi, . . . ,Vd are affine functions, 
i.e. we have 

Vi-.S-+R n , x^HiX + hi, 
where Hi G R" xn and hi G M. n . A process X which solves the stochastic differential 
equation of type 

d 

dX t = Vi{X t _)dL t>i , X = xeS, (4.3) 
i=i 

in S and which leaves S invariant is m-polynomial */ijif||>i < 00 ■ V 

V\, . . . , Vd are linear, then we even have A(Pol k (S)) C Polfc(S') for all < k < m 
(in contrast to A(Pol< k (S)) C Pol< fe (S%). 

Proof. For u G C 2 (S) C T>(A) and general Lipschitz continuous functions V\, . . . , Vd 
the infinitesimal generator of X is given by 



1 - 

Au{x) = - ((W ■ ■ ■ V d {x))a{Vx{x) . . . V d {x))') kl 



k,l=l 
d 



d 2 u(x) 
dx k dxi 



(^Vi{x)h,Vu{x)j 



u(x) 



53^(a:)ift,Vtt(ar)^ n{dy). 



This corresponds again to the situat ion of Condition B w ith g{x,y) = H(y)x 
Kv) — Z)i=i HiViZ + hiVi- Theorem 



3.6 



and Remark 

slight generalization of (4.3 1, then yield the assertion. 
If V x , 



3.9 



which actually treats a 



, Vd are linear, the entries of the diffusion matrix are quadratic and the 
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drift is linear. This compensates the decreasing degree of the second and the first 
derivatives respectively. In the jump part we only get terms of order m. Thus, in 
this case A maps Pol m (5') into itself. □ 

Example 4.5 (Jacobi process). Another example of a polynomial process is the 
Jacobi process (see |10] ) which is the solution of the stochastic differential equation 



dX t = -[3{X t - 6)dt + ay/Xt(l-X t )dB t , X = xe [0, 1], 

on S = [0, 1], where 9 £ [0, 1] and /3, a > 0. Its generator is a self-adjoint operator 
satisfying the following eigenvalue equation 

AQ n (x) = X n Q n (x), neN, 

where the eigenfunctions Q n are the Jacobi polynomials and the eigenvalues X n are 
given by ~^-n(n — 1 + |f )■ 

This example can be extended by adding jumps, where the jump times correspond 
to those of a Poisson process N with intensity X and the jump size is a function of 
the process level. Indeed, if a jump occurs, then the process is reflected at ~ so that 
it remains in the interval [0,1], i.e. we have 

dX t = -p(X t - 9)dt + ay/X t (l-X t )dB t + (1 - 2X t )dN u X = x e [0, 1], 

whose generator is given by 

, 1 2/ /, ..d 2 u(x) n .du(x) , , ,„ s , ss 

Au = -a 2 (x(l - x)) , \ ' - j3ix - 9)—^- + X(u(l - x) - u(x)). 
2 dx z dx 

In terms of Condition B, we have here g(x,y) = —2yx + y and fJ-(dy) = XS\(dy). 



Example 4.6 (Pearson diffusions). The above example 4-5 (without jumps) as well 
as Ornstein-Uhlenbeck and Cox-Ingersoll-Ross processes, all of them with mean- 
reverting drift, can be subsumed under the class of so called Pearson diffusions 
which are the solutions to SDEs of the form 



dX t = -f3(X t - 0)dt + ^(a n X? + a w X t + a)dB u X = x, 

where (3 > and an, aio and a are specified such as the square root is well defined. 
Forman and S0rensen [5] give a complete classification of the different types of the 
Pearson diffusion in terms of their invariant distributions. 



5. Applications 

By Theorem |2.3| we know that there exists a linear map A such that moments 
of m-polynomial processes can simply be calculated by computing e tA . Indeed, by 
choosing a basis (ei, . . . , ejy) of Pol< m (S') the matrix corresponding to this linear 
map which we also denote by A — ,j=i,...,jv can be obtained through 

N 

Aei — Aij ej . 

3=1 

Writing / as X)fc=i a kS-k, we then have 

(N \ N 

a k e k = Y a ke tA e k = (ai, . . . , a N )e tA {e x , e N )', (5.1) 
k=l ) k=l 



POLYNOMIAL PROCESSES AND APPLICATIONS TO MATHEMATICAL FINANCE 17 



which means that moments of polynomial processes can be evaluated simply by 
computing matrix exponentials. 

By means of the one-dimensional Cox-Ingersoll-Ross process 

dX t = (b + pX t )dt + (TyfXtdWt, b,aeR+, (3 e R, 

we exemplify how moments of order m can be calculated. The generator is given 
by 

Au{x) 



n d?u(x) ., n .du(x) 



Applying A to [x , x x m ) yields the following (m + 1) x (m + 1) matrix 



.4 



/ ... 

b (3 ... 

2b + a 2 2(3 

36+3cr 2 3/3 



V o .. 

Hence, E[(Xf) fe ] = P t x k 



(0, 



1. 



nib - 
,0)e tA (a; ,. 



m(m — 1) 2 
2 G 



m(3 J 

x m y. 



Remark 5.1. Note that A is a lower triangular matrix, whose eigenvalues are the 
diagonal elements. Since in this case they are all distinct, the matrix is diagonal- 
izable. This holds true for all one- dimensional affine processes since the (k + l) th 
diagonal element is given by 

where the integral part can only appear if the process is supported on M + . Recall 
the notation hiq from Condition A. Of course there are many efficient algorithms 
to evaluate such matrix exponentials. 

Remark 5.2. If n > 1 one has to apply well-known techniques from linear algebra 
(such as in FEM) in order to enumerate efficiently basis o/Pol< m (5) and to exploit 
sparsity properties of A. 

5.1. Moment calculation - Generalized Method of Moments (GMM). In 

view of this easy and fast technique of moment calculation for polynomial processes, 
the Generalized Method of Moments (GMM) qualifies for parameter estimation and 
thus, for model calibration. The implementation of a typical moment condition of 
the type 



/ X^X^-EiX^X™^ 



f(X t ,9) = 



n^rrii € N, l<i<q, 



where 9 is the set of parameters to be estimated, is simple since 

E[X?X%_ a ] =E[X?E[X%. a \X t ]] 

can also be computed easily. In the case of one-dimensional jump-diffusions, 
Zhou already uses this method for GMM estimation. 
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5.2. Pricing - Variance reduction. The fact that moments of polynomial pro- 
cesses are analytically known also gives rise to new and efficient techniques for 
pricing and hedging issues. 

Let X be an m-polynomial process and G : S — > R" a deterministic bi-measurable 
map such that the price processes are given through 

S t = G(X t ) 

under a martingale measure. Typically G = exp, if X are log-prices. We denote by 
F = 4>(St) a bounded measurable European claim for some maturity T > 0, whose 
price at t > is given by the risk neutral evaluation formula 

pf =E[<l>(S T )\ft] = E[(<t>oG)(X T )\F t ]. 

Obviously, claims of the form 

F = foG- 1 (S T ) (5.2) 

for / <E Pol< m (5) are analytically tractable, as we have 

pf = E[(fo G-^St)^] = Pr-tfiG-^St)) = e^ A f{G-\S t )) 

for < t < T, where A is the previously defined linear operator on Pol< m (S f ). 
The sensitivities of the price process with respect to the factors of X can then be 
calculated by 

Vpf = VPr-tfiG-^St^VG-^St). (5.3) 



Although claims are in practice not of form (5.2), the explicit knowledge of the 
price of polynomial claims can be used for variance reduction techniques based on 
control variates. Instead of using the estimator 



^=^I>°G)(A^) 

i=l 

in a Monte-Carlo simulation, where X^, ■ ■ ■ ,X^, are L samples of Xj>, we can use 



yJ2 ((* G)( x t) - (/TO ~ E[f(X T )])) 



L 



where / G Pol< m (S l ) is an approximation of <f> o G and serves as control variate. 
Both estimators are unbiased and the second clearly outperforms the first since 
Var (fto) < Var (^o), where the ratio of the variances depends on the accuracy of 
the polynomial approximation. 



It is worth mentioning that the previous pricing algorithm has also important 
consequences for hedging, since the greeks for "polynomial claims" F — f(Xx) can 
be explicitly and efficiently calculated, again by matrix exponentials: the matrix 
exponential calculates the coefficients of the polynomial x i— > E[f(Xx)], calculating 
a derivative of this polynomial is then an algebraic question. These considerations 
are then applied in equation ( |5.3[ ). Now assuming a complete market situation for 
the real-world claim 4>(St) — 4>° G(Xt), i.e. there is a trading strategy rj such that 

<t>(S T ) = E[<p(S T )}+ I rjfdSt, 
Jo 
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we can conclude that 

cf>(S T ) ~ f(X T ) = EMSt)] ~ E[f(X T )] + - Vpf ) • dS t . 

Jo 

Therefore, if we assume that </>(Sr) — f(Xr) has a small variance, then also the 
stochastic integral representing the difference of the cumulative gains and losses of 
the two hedging portfolios, namely the one built by the unknown strategy rj and 
the one built by the known strategy Vpf , is small. 

Of course, this method of variance reduction can be applied to all polynomial 
processes, including all kinds of examples mentioned in section [I] However, afhne 
models for which the generalized Riccati ODEs (see [5]) cannot be explicitly solved 
are of particular interest. This is illustrated in the example below. 

Example 5.3. The following affine stochastic volatility model comprising two vola- 
tility factors is a modification of a model initially proposed by Bates [2] . The price 
process is specified as St — Soe Xt with dynamics 

,{*t\ /r-a-^-A 1 C/ t /(e«i-l)F 1 (d€)-A 2 V t /(e«i-l)j' 2 (d|))\ 

d I u t = h-PuUt+PizVt dt 

\ Vt / V b 2 ~f3 2 2Vt+f32iUt J 

I sJUl JV t \ /dB tA \ 

+ viPiVlh viy/l-pjVTTt ° dBtl 

\ o- 2 p 2 VVi <t 2 t/1-P%V%J \dBt'^J 

+ dZ t .i + dZ tj2 , 

where for k — 1,2, Z^ are pure jump processes inlx M. 2 ^ with linear jump arrival 
intensity \\U, X2V and trivariate jump size distribution F^. The usual way to obtain 
the price of European options is to solve the following Riccati equations 

d(p{t,x,u,v) , 

■ = rx + Oiipi(t, x, u, v) + 02ip2{t, x, u, v), 



dt 

II 1U 1 j 

(x - x) - Pii%l>i{t,X,U,v) + 0211p2{t,X,U,V) 



dlp\(t,X,U,v) 1/2 



dt 2 



+ \ 2 G \i>\ l K t i X : M : V ) + Pl<J\Xi>\(t, X, U, V) 



+ A i |^ y ^ ^ e ^i+Vi(t,!<:,«,f)6+'A2(t,s,«,u)?3 _ Fi(d£) 
-x f (eti-l)Fi(#) V 

JrxR% J 

dlp2{t,X,v) 1/9 s n 1 \ nil \ 

= -(X - X) - P22W2(t,X,U,V) + Pi2lpl(t,X,U,V) 



dt 2* 



+ q ^^*! X ' U > + P2^2X1p2{t, X, U, v) 



+ A 2 |^ J ^ ^ e a:5l+i/'l(*:a:>«,i')?2+i/ , 2(t,a;,M,f)C3 _ ^ i^^O 

-x f (efr-ljFsCdO). 
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and to apply Fourier pricing methods as suggested in |6J. As in this case explicit 
solutions of the Riccati equations are not available, our approach to use an approx- 
imating polynomial as control variate in the Monte-Carlo simulation is particularly 
expedient. For the calculation of the matrix exponential yielding the price of the 
polynomial claim, we only need to apply the infinitesimal generator 

8f , ,u o .. , o „.\ df , , t , fl df 

dv 

'u v\d 2 f ua\d 2 f va 2 d 2 f d 2 f d 2 f 



(U v\ Of u T 

T ~ 2 ~ 2 / 8 + (61 ~ /9uU + /9l2W - ) a + (&2 ~ ^ + /32lU) 
/it u\ 
+ (2 + 2) 



9a; 2 ' 2 <9u 2 ' 2 <9u 2 ' <Tl ^ 1 9u<9a; 2 ^ 2V ~dvdx 

6 



+ A lM y ^/( a; + ^, u + ^, u + ^)_/( a; ,u ) t;)-(e ei -l)^jF 1 (dC) 

+ A 2W y ^/(x + Ci,w + 6^ + 6)-/(^^w)-(e ei -l)^F 2 ( f iO 

polynomials of the form x kl u k2 v ki , ki + k 2 + k 3 < m in order to determine the 
matrix A. The price of the approximating polynomial claim is then simply calculated 



by means of (5.1 1 



In order to present the positive impact of our variance reduction method graph- 
ically, we implemented the following simplified version of the above model (indeed 
the original Bates model ) with dynamics of the form 

V t J ~ \ b-0V t ) at 

W t \ / dB tA \,(dZ t 

apVV t o-^T^VVt J \ dBt,2 J V 
where Z is a pure jump process in K with jump intensity Xv and jump size distribu- 



tion F. Figure 5.2 illustrates the comparison between the Monte Carlo simulation 
for European call prices with and without variance reduction. 

Our variance reductions technique can also be successfully applied to derivatives 
involving several assets, provided that their dynamics are described by a polynomial 
process. This can be done simply by approximating European payoff functions 
depending on several variables with multivariate polynomials. 
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